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Context. Ground-state OH maser emission from late-type stars is usually polarized and remains a powerful probe of the magnetic 
field structure in the outer regions of circumstellar envelopes if observed with high angular and spectral resolutions. Observations in 
all four Stokes parameters are quite sparse and this is the most thorough, systematic study published to date. 

Aims. We aim to determine polarization properties of OH masers in an extensive sample of stars that show copious mass loss and 
search for candidate objects that are well-suited for high angular resolution studies. 

Methods. Full-polarization observations of the OH 1612 and 1667 MHz maser transitions were carried out for a sample of 1 17 AGB 
and post-AGB stars. Several targets were also observed in the 1665 MHz line. 

Results. A complete set of full-polarization spectra together with the basic polarization parameters are presented. Polarized features 
occur in more than 75% of the sources in the complete sample and there is no intrinsic difference in the occurrence of polarized 
emission between the three classes of objects of different infrared characteristics. The highest fractional polarization occurs for the 
post-AGB+PN and the Mira+SR classes at 1612 and 1667 MHz, respectively. Differences in the fractional polarization between the 
sources at different evolutionary stages appear to be related to depolarization caused by blending. The alignment of the polarization 
angles at the extreme sides of the shell implies a regular structure of the magnetic field of a strength of 0.3-2.3 mG. 
Conclusions. Polarized OH maser features are widespread in AGB and post-AGB stars. The relationship between the circular and 
linear fractional polarizations for a representative sample are consistent with the standard models of polarization for the Zeeman 
splitting higher than the Doppler line width, whereas the polarized features are the a components. 
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1. Introduction 

The ground-state OH maser lines at 1612, 1665 and 1667 MHz 
are frequently detected in late-type stars that show mass loss at 
high rates. A recent compilation by Engels & Bunzel (2008) 
gives 2297 circumstellar masers of which 1945 show the clas- 
sical double-peak profile at 1612 MHz. The polarization infor- 
mation is not available or neglected in the interferometric sur- 
veys of the Galactic plane (Sevenster et al. |1997a|[T997Hl200TT l 
in numerous single-dish studies of IRAS colour-selected objects 
(e.g. Eder et al. [l988l Le Squeren et al. lTWTl te Li ntel H ekkert 
et al. l 199 It and in earlier blind surveys (e.g. Bowers 1978; Baud 
et al. |1979al |1979bll. E arlier studies (e.g. Wilson et al. [T970l 
Wilson & Barrett l 1972b suggested that the 1612 MHz circum- 
stellar emission is essentially unpolarized but the observations 
with a spectral resolution better than O.lkms -1 revealed sev- 
eral circularly polarized features (Cohen et al. l 19871 Zell & 
Fix [T99"n >. The OH emission in the 1665 and 1667 MHz lines 
shows a substantial degree of circular polarization up to 50-80% 
(Wilson et al. l 19721 Claussen & Fix 1982 ) but linear polarization 
occurs in only ~10% of sources (Olnon et al. |1980t . Highly po- 
larized OH maser features are observed in flaring Miras (Reid 
et al. [T977l Fix [l979l Eto ka & L e Squeren [T997] l. hypergiant 
IRC+ 10420 (Benson et al.[T979l and post-AGB star OH17.7- 
2.0 (Szymczak & Gerard2005). There is a general consensus 
that the polarization measurements of OH maser emission can 
be properly explained with the classic model of Goldreich et al. 



* Appendices (Table A. 1) and (Figure B.l) are only available in the 
electronic edition of the journal at http://www.aanda.org 



(1 1973b in which the magnetic field plays a key role (e.g. Gray & 
Field [T996l Elitzur [l996l 

The presence of circular polarization is itself evidence for 
Zeeman splitting (Deguchi & Watson 1986 ) and, despite the dif- 
ficulties in identifying classical Zeeman patterns, provides a di- 
rect estimate of the strength and direction of the magnetic field. 
From observations of the linear polarization we can draw conclu- 
sions about the orientation of the circumstellar magnetic field. 

For some OH stellar sources, high-sensitivity observations 
with ~0.1kms~' spectral resolution have shown an amount 
of circularly and linearly polarized features large enough to 
measure the strength of circumstellar magnetic fields and to 
map their structure with high accuracy using the interfero- 
metric technique (e.g. Chapman & Cohen 1986; Szymczak et 
al. [T998l 120011 Bains et al. l2003l Amiri et al. l2010l ). These case 
studies revealed the presence of highly ordered magnetic fields 
of a strength of a few mG at distances of a few thousands stel- 
lar radii. Full -polarization measurements of OH circumstellar 
masers with high sensitivity and spectral resolution are scarce. 
A small sample of Miras (Claussen & Fix 1982) and a mod- 
erate sample of proto-planetary nebula candidates (Szymczak 
& Gerard 2004) only were studied up to date and no compre- 
hensive picture of the polarization properties of circumstellar 
OH masers was obtained. We hereby extend those studies for a 
sample containing oxygen-rich objects that evenly populate the 
whole evolutionary sequence in the IRAS two-colour diagram 
(van der Veen & Habing 1988 ). There is evidence that the IRAS 
colours increase through the sequence from the semi-regular 
variables (SRs), long period variables of Mira types, OH/IR 
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Fig. 1. IRAS colour-colour diagram for the sample studied. The 
regions defined by van der Veen & Habing ( 1988 ) for different 
types of objects are marked by dotted lines. The three classes of 
objects discussed in the text can be distinguished by the symbols. 



cerned here only with these three groups of objects. Note that 
the latter group was included in a sample studied by Szymczak 
& Gerard (2004) but the spectra used there were generally taken 
at other epochs. One hundred and five out of 117 targets have 
good flux quality at the 12, 25 and 60/im IRAS bands, and we 
show their two-colour diagram in Figure Q] The colour is de- 
fined as [A/] - [Aj] = -2.5log{S(Ai)/S{Aj)), where S(^ ; ) is the 
uncorrected flux in the IRAS /I, band. The colours of OH/IR and 
post-AGB objects overlap, partly because of inaccurate classifi- 
cation criteria of the transition objects between AGB and post- 
AGB phases. Figure Q] shows that our targets sample all types of 
oxygen-rich stars with a circumstellar envelope associated with 
OH maser emission. The sample is therefore appropriate to prop- 
erly address whether the polarization parameters of OH masers 
depend on the characteristics of the central stars/envelopes or are 
dominated by propagation effects. 

A comparison of the 1612 MHz peak flux density for 105 
detected sources in the sample with that for 808 objects in the te 
Lintel Hekkert et al. ( 119911 ) catalogue indicates that our sample 
probes the brightest sources; the median values are 8.32 and 0.86 
Jy, respectively (Table [1]). We note, however, that the older ob- 
servations may underestimate the peak flux because of frequency 
dilution. 



objects, proto-planetary nebulae (PPN) and planetary nebulae 
(PN). The colours probably reflect the evolution of the circum- 
stellar shell from optically thin to optically thick and the location 
of an object in the diagram could be also related to the age, initial 
mass and possibly metallicity (e.g. Bedijn 1987; van der Veen & 
Habing 1988 1 Likkel 1989). The observations of a representative 
sample reported here attempt to examine whether the OH maser 
polarization depends on the properties of the central star and its 
envelope or whether it is caused by propagation effects. 

2. Data 

2.1. The sample 

The initial source list consisted of 152 late-type stars at 6 <; -39° 
of which 147 objects were taken from Benson et al.'s catalogue 
d 19901 ). the rest from te Lintel Hekkert et al.'s survey (11991I ). 
Several sources were discarded after the initial test observa- 
tions because their total peak fluxes were lower than 0.15 Jy at 
both frequencies. The final sample consisted of 117 stellar OH 
sources that clearly exhibit either the 1612 MHz satellite line 
or/and the 1665/1667 MHz main lines. Targets with different in- 
frared characteristics were chosen to probe the whole evolution- 
ary sequence of oxygen-rich objects in the IRAS two-colour di- 
agram (van der Veen & Habing 1988 ). The location of an object 
in this diagram depends on the circumstellar envelope thickness, 
age and initial mass. It is assumed that the visible semi-regular 
and Mira variables with relatively thin envelopes evolve towards 
OH/IR objects with thick envelopes, and when their pulsations 
cease, they leave the tip of the asymptotic giant branch and be- 
come post-AGB objects and eventually young planetary nebu- 
lae (Habing 1996 for review). Alternatively, OH/IR stars have 
higher progenitor masses than Miras (Likkel 1989) and may 
come from an intrinsically more massive population (Whitelock 
et al. 119941 ). These categories of objects differ in pulsational 
and mass-loss properties and can be accompanied by OH maser 
emission. 

There are 36 SRs and Miras, 57 OH/IR objects and 24 post- 
AGB and PN candidates in the sample (Table A.l). We are con- 



2.2. The observations 

The data were obtained with the Nangay Radio Telescope (NRT) 
during several runs between February 2002 and December 2008. 
This transit instrument is the equivalent of a dish of 93 m in di- 
ameter and the half-power beam-width is 3^5(RA)xl9'(Dec) at 
1.6 GHz. At Dec = 0°, the beam efficiency is 0.65 and the point- 
source efficiency is 1.4KJy~', the system temperature is about 
35 K. More technical details on the radio telescope and its up- 
grade are given in van Driel et al. ( 1996|. 

Initially, we observed 152 sources in the 1612.231 and 
1667.359MHz OH lines, 38 of which were also observed at 
1665.402 MHz (including 31 sources that were not detected at 
1612.231MHz). We used a correlator configuration providing 
eight banks of 1024 channels. To cover all OH maser profiles 
easily, we used two spectral resolutions, ~0.14kms _1 for about 
two thirds of the targets and ~0.07kms _1 for the rest. 

Two transitions were simultaneously observed and we reg- 
istered the two orthogonal linear polarizations and two oppo- 
site circular polarizations for each transition. The system di- 
rectly provided three of the four Stokes parameters, namely 
I, Q and V, while the fourth parameter U was extracted by a 
horn rotation of 45°. The linearly polarized flux density, p = 

V2 2 + U 1 , fractional linear polarization, «l = p/I, fractional 
circular polarization, mc — V/I and polarization position angle, 
X = 0.5xtan~ l (U/Q), were derived from the Stokes parameters. 
The methods of observation, polarization calibrations, data re- 
duction and estimation of uncertainties on the Stokes parameters 
were the same as described in Szymczak & Gerard (2004, Sect. 
2). Briefly, we observed in frequency switching mode, using a 
noise diode to measure the gain of each polarization bank at the 
start of each scan. The absolute and relative accuracies of gain 
were ~5% and ~1%, respectively. The instrumental polarization 
was regularly checked with observations of W12 and W30H and 
the absolute flux density scale was accurate to within 7-8%. The 
error in the polarized intensity caused by the polarization leak- 
age between the orthogonal feeds was about 2%. The typical 
integration time for each horn position was about 20 min. This 
resulted in a rms noise level in the Stokes / of about 35 mJy for 
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Table 1. Comparison of mean and median peak flux densities in 
our sample and that of te Lintel Hekkert et al. dl9911 l. 



sample number frequency mean(SD a ) median 

of sources (MHz) (Jy) (Jy) 



this paper 105 1612 31.70 (8.12) 8.32 

85 1667 3.59 (0.74) 1.21 

te Lintel Hekkert 808 1612 4.37 (0.84) 0.86 



3.2. Occurrence of polarized emission 

The detection of polarization is limited by the noise level of 
0. 15 Jy in any Stokes parameter and the minimum believable 
fractional polarization of 2%. Thus, our sample is complete for 
objects that are higher than 2% polarized and whose total in- 
tensity exceeds 7.5 Jy. The weaker sources only have detectable 
polarization if they are more strongly polarized. Table [3] summa- 
rizes the detection statistics for elliptical, circular and linear po- 
larizations towards the sources divided into three groups depend- 
ing on the occurrence of the 1612 MHz or/and 1667 MHz lines. 
Column 2 gives the number of detected sources and columns 3-5 
list the number of sources with features of given types of polar- 
ization. The emission at 1612 and 1667 MHz was detected in 105 
and 85 targets, respectively. Polarized features were detected in 
75% (79/105) and 46% (39/85) of the sources at 1612MHz and 
1667 MHz, respectively. This is because the 1612 MHz emission 
is generally stronger than the 1667 MHz emission. The study in- 
dicates that elliptical polarization dominates at 1 6 1 2 MHz. There 
is also a high detection rate of linearly polarized emission (58% 
at 1612 MHz and 26% at 1667 MHz) in the sample, never re- 
ported previously (Cohen 1989, for review). Surprisingly, we de- 
tected 14 sources with linearly polarized features without appar- 
ent circular polarization (see Sect.3.5). 

3.3. Fractional polarization 

The polarization properties of maser spectra are commonly de- 
termined relative to the peak flux density of Stokes / parameter. 
However, we found that the velocities of the Stokes / peak do 
not coincide with those of the Stokes V and p peaks for 62% and 
65% of the sources, respectively. 

The median velocity differences between the peaks of Stokes 
/ and Stokes V and p are 0.14 and 0.27 km s~', respectively. 
Generally, the most strongly polarized features coincide with the 
Stokes / peaks within 2-3 spectral channels. However, for a few 
sources these velocity differences are as large as 5.2kms _I in 
IRAS 191 14+0002, Ukms 1 in OH12.8-1.9 and 0.9 km s" 1 in 
RR Aql, then the V and p peaks are distinct features. Velocity 
separations could be caused by Zeeman splitting (see Sect. 3.6). 
Blending of differently polarized components along the line of 
sight would also cause depolarization relative to the total inten- 
sity. In the following sections, we perform separate statistical 
analyses for the fractional polarization at the Stokes /, Stokes 
V and p peaks. We compare the results for the three classes of 
objects defined in Table A.l. 

3.3.1 . Measurements relative to the Stokes / peak 

Figure|2]summarizes the polarization properties of OH emission 
for the three classes of objects in the sample, and Table |4]gives 
the mean and median values of circular and linear fractional po- 
larization. The fractional circular polarization in the OH sources 
is commonly lower than 20% and is as high as 60-100% in a 
few objects. In contrast, the fractional linear polarization is usu- 
ally lower than 10% and reaches a slightly higher value in a few 
sources. The differences in the fractional polarization distribu- 
tions between the three classes and the two lines were investi- 
gated using a Kolmogorov-Smirnov test. In the post-AGB+PN 
objects the mean value of mc at 1612 MHz is higher (p < 0.01) 
than that for the other two classes. The mean degree of circu- 
lar polarization at 1667 MHz of the Mira+SR objects is higher 
(p < 0.01) than that of the OH/IR and post-AGB+PN objects. 
There is no significant difference in the distribution of be- 



Notes. " standard deviation 

Table 3. Number of sources with polarized features at the OH 
lines. 



Line 


Number of 


Number of sources with features 


(MHz) 


detection 


elliptical 


circular 


linear 


1612 only 


32 


10 


4 


3 


1667 only 


12 


3 


7 





1612/7667 


73 


52 


16 


7 


1667/7672 


73 


15 


8 


4 


total 


1612 


105 


58 


20 


10 


1667 


85 


18 


15 


4 



0.14 km s~' spectral resolution. The median 3cr rms noise levels 
of the Stokes Q and U spectra were 0.15Jy. In those features 
that showed polarization, the polarization was only considered 
believable if either or both of the linearly or circularly polarized 
fluxes were higher than 0. 15 Jy. 

The data for our sample have been collected over nearly 
seven years. Each target was observed at least three times on an 
irregular basis but there are several sources that were observed 
over a 3-7 years period so that the light curves could be easily 
obtained. Here we only report data for the epoch during which 
the maximum flux density occurred. Results of the whole long- 
term observations will be published in a separate paper. 



3. Results and analysis 

Table A.l contains a list of the 117 detected sources and their 
polarization parameters. Full-polarization spectra are shown in 
Figure B.l. 

3.1. Effect of spectral resolution 

We tested whether the fractional polarization decreases as the 
spectral resolution is degraded by observing a few targets at 
both spectral resolutions. For the two sources OH17. 7-2.0 and 
OH138.0+7.2 the maximum values of mc and at 1612 and 
1667MHz decreased by less than 0.2% and 0.5%, respectively, 
when the resolution was lowered from ~0.07 to ~0.14kms _1 . 
Additionally, the 1612 and 1667MHz spectra of ~0.07kms _1 
resolution with detected polarized features were smoothed to a 
resolution of ~0.14kms _1 . The average and median values of 
the fractional polarizations for the strongest features before and 
after the smoothing are given in Table [2] No significant differ- 
ence in the fractional polarizations is found for both observed 
lines. These two tests clearly indicate that our data taken with 
the two different spectral resolutions can be treated together in 
the subsequent analysis. 
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Table 2. Average and median values of the fractional circular (mc) and linear (mi_) polarization for the strongest features at the two 
spectral resolutions. 



Line 


Fractional 


Number 


Mean(SD)° 


Median" 


Mean(SD) 1 


Median 1 


(MHz) 


polzn. 


of sources 


(%) 


(%) 


(%) 


(%) 


1612 


\mc\ 


11 


7.5(1.8) 


5.8 


8.2(2.3) 


5.5 




m L 


12 


4.5(0.6) 


4.0 


4.2(0.6) 


3.9 


1667 


\mc\ 


12 


27.0(9.4) 


10.1 


27.0(9.0) 


10.2 




m L 


6 


7.1(2.6) 


4.2 


7.5(2.6) 


5.6 



Spectral resolution: ~0.07 km s 
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Fig. 2. Polarization properties of OH emission in three classes of late-type stars. The flux density measurements are made at the 
Stokes I peak. Top: Histograms of circular (mc) and linear (tol) fractional polarizations for each class of objects at both OH lines. 
Middle: Mean and median values of the circular and linear fractional polarizations for three classes of objects. The number of 
objects for each class and OH line is given. Bottom : Circularly polarized peak flux density (Stokes | V|) (left) and linearly polarized 
peak flux density (parameter p)(right) versus the total peak flux density (Stokes /) for the two OH lines. The coloured (online) 
symbols represent classes of objects while the grey symbols mark the sources with polarized emission lower than the individual 
detection level (3<r). The dotted lines show the mean upper limit 0.15 Jy. Some objects have much lower individual upper limits 
because they were observed longer. The dashed lines labelled 2, 10, 30 and 100 indicate the percentage of polarization. 

Table 4. Mean and median values of fractional polarization determined at the velocity of the Stokes / peak flux density. 



Class 


Line 


number 


\m c \(%) 


number 


m L (%) 




(MHz) 


of sources 


mean (SD) median 


of sources 


mean (SD) median 



Mira+SR 


1612 


9 


8.70(1.14) 


4.62 


9 


3.89(0.66) 


3.89 




1667 


14 


24.98(4.95) 


10.13 


6 


7.70(0.66) 


5.14 


OH/IR 


1612 


33 


5.81(0.48) 


4.28 


25 


4.70(0.29) 


4.34 




1667 


11 


8.28(1.36) 


4.55 


8 


5.71(0.59) 


5.05 


post-AGB+PN 


1612 


10 


24.71(1.17) 


14.14 


11 


5.50(0.44) 


4.52 




1667 


3 


11.30(0.82) 


3.87 


2 


5.60(0.24) 


5.60 
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Fig. 3. Comparison of the mean fractional polarizations obtained 
with two different methods (Tables|4]and|5j- The mean values of 
niL and mc correspond to the rectangle centres and the standard 
deviations are marked by the borders. The origins of the arrows 
mark the mean fractional polarizations estimated at the Stokes / 
peaks, whereas the arrowheads mark the mean fractional polar- 
izations estimated at the Stokes V and p peaks. 
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Fig. 4. Histograms of all sources versus the integrated Stokes 
/ flux density for weak and more strongly polarized 1612 and 
1667 MHz lines. The fractional polarizations Mc and Ml are the 
ratios of the integrated flux densities. Note that there are statis- 
tically significant differences (p < 0.01) between the two distri- 
butions in each panel. 



tween the three classes of objects. There are 10 Mira+SR stars, 
21 OH/IR objects and 8 post-AGB+PN sources with I p >10Jy 
without detectable polarization (Fig. |2j. 



3.3.2. Measurements relative to the Stokes V and p 

The velocity offset of the Stokes V and p peaks from the Stokes / 
peaks, as noted in Sect. 3.3, is most pronounced for the Mira+SR 
stars where it is observed in 70-75% of the objects. The result of 
the analysis of the fractional polarization based on the measure- 
ments concerning the peak fluxes in the V and p parameters is 
summarized in Table [5] Generally, the mean and median values 
of mc and are 2-3 times higher than the values measured 
at the Stokes / peak (Table |4). The increase in mc is greater at 
1667 than at 1612 MHz. There is a remarkable increase of both 
\m c \ and m L for the 1667 MHz Mira+SR objects (Fig. [3>. The 
trends described in Sect. 3.3.1 remain valid. We also measured 
the dispersion of mc and tol in consecutive spectral channels 
with polarized signal. It appears that at 1612 MHz the dispersion 
of mc for the Miras+SR and post-AGB+PN stars is nearly twice 
higher than in the OH/IR objects. The Mira+SR variables also 
show the highest dispersion of mc at 1667 MHz. The present 
analysis clearly shows that for the majority of sources the most 
highly polarized features are weaker than the line peaks. In the 
following analysis of the polarization properties, the measure- 
ments relative to the Stokes V and p are used because they are 
the best estimates of the maximum values of the polarization pa- 
rameters in the sources (Tables [4] and [5]). 



3.3.3. Ratios of the integrated flux densities 

There are several sources in the sample that exhibit complex pro- 
files in the Stokes V and p so that their polarization properties 
can only be properly determined using the integrated flux densi- 
ties of the Stokes parameters. The number of sources versus the 
integrated Stokes / flux density for both OH lines are shown in 
Fig.il Here M c = VintMnt and M L = pintAnt, where / int , V int and 
Pint are the integrated Stokes /, V and p flux densities, respec- 
tively. A distinction is made between weakly (Mc, Ml <3%) and 
more strongly (>3%) polarized sources. There is evidence for 
higher circular and linear polarization in the 1667 MHz sources 
with /j nt < 1 0-20 Jy km s _1 . Less polarized 1667 MHz sources are 
usually the brightest in the sample. For the 1612MHz line there 
is a slight excess of weak (<10 Jy kms~') sources with Mc >3%. 
For the majority of the 1612 MHz sources Mc and Ml do not de- 
pend on the integrated Stokes / flux density. Here the defined ra- 
tios estimate the net circular and linear polarizations. The net 
polarization would be zero for a uniformly filled maser shell 
with symmetric magnetic and velocity fields and density gra- 
dient. Hence, the trend to higher integrated fractional circular 
polarization at lower integrated total intensity may suggest that 
these objects have less homogeneous velocity and/or density dis- 
tributions. 

3.4. Highly polarized sources 

Fractional polarizations referred to in this section as mc and m^ 
are determined at the velocities of the V and p peaks. Table [6] 
lists the sources with fractional circular and linear polarizations 
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Table 5. Mean and median values of fractional circular and linear polarization determined at the velocities of the Stokes V and p 
peak flux densities, respectively. 



Class 


Line 


number 


|m c l(%) 


number 


m L (%) 




(MHz) 


of sources 


mean (SD) median 


of sources 


mean (SD) median 



Mira+SR 


1612 


14 


12.62(1.68) 


5.57 


12 


7.13(1.15) 


4.07 




1667 


18 


39.31(8.76) 


26.53 


8 


11.34(0.93) 


8.56 


OH/IR 


1612 


38 


7.37(0.68) 


6.13 


37 


5.46(0.40) 


4.75 




1667 


12 


11.54(1.71) 


8.70 


9 


11.32(1.41) 


8.77 


post-AGB +PN 


1612 


18 


23.41(3.47) 


13.47 


11 


6.79(0.44) 


5.19 




1667 


6 


24.51(3.42) 


26.28 


4 


4.23(0.24) 


2.96 



Table 6. Sources with highly polarized OH features. The frac- 
tional polarizations for the Stokes V and p peaks, V p and p p are 
given. For the Stokes / peak the fractional polarization is also 
given when we were able to reliably determine it. 



Name Line m c (V p ) m c (I v ) m L (p p ) m L (I p ) 

(MHz) (%) 



greater than 20% and/or 10%, respectively. There are 26 sources 
with mc >20% and 17 sources with ;«l >10%. 

This subsample contains 10 and 21 objects with the 1612 and 
1667 MHz lines, respectively. About two thirds of the 1667 MHz 



highly polarized sources are the Mira+SR variables. Fig. 
shows the relation between the circular and linear fractional po- 
larizations for the whole sample. The majority of moderately 
polarized sources, i.e. with well-determined \mc\ <25% and 
mL <10% lies below the line of form \mc\ = 2.1 + 1.5. 
There are a few groups of outliers; the group of Mira+SR 
objects with strongly (\mc\ >75%) polarized 1667MHz fea- 
tures without measurable linear polarization, the two OH/IR 
objects (OH17.7-2.0 and IRAS 16342-38 14) with features of 
\m c \ >80% and m L £15%, the two Miras (GY Aql and UXCyg) 
with features of high linear and circular polarization, the group 
mainly of OH/IR objects with features of >10% and low or 
marginal circular polarization. We notice that the sources with 
highly polarized features are usually the Mira+SR objects and 
a few OH/IR objects showing outburst activity in OH maser 
lines (e.g. UXCyg - Etoka & Le Squeren d20D0l OH17.7-2.0 
- Szymczak & Gerard d20"05V ). 

3.5. Polarization angle arrangement 

There are 29 1612 MHz sources with linearly polarized emission 
in both parts of the double peaked profile. Several exemplary p 
spectra with the superimposed position angle of linear polariza- 
tion are shown in Fig. [6] Most of the sources show only weak 
(<15°) variations of the polarization angle from channel to chan- 
nel. Table [7] gives the average polarization angles for the four 
most blue- and red-shifted channels in the 1612 MHz p profiles 
(Xb, Xr) ar, d me difference (A^br)- We found that A^br is lower 
than 20°for two thirds (20/29) of the sources and is only higher 
than 70° in two sources. For the whole subsample the absolute 
average and median values of A^br are 19°3+4.0 and 1 1°6, re- 
spectively. For a few sources observed at several epochs span- 
ning 6-8 years, no systematic variations of A^bs higher than 10° 
were found. 

In the standard circumstellar shell model (Reid et al. 1977 ) 
the 1612 MHz double peaks come from two very compact re- 
gions of opposite sides of the shell. VLBI observations of 
the archetypal OH/IR object OH127. 8+0.0 proved that the 
1612 MHz emission of ~0.5kms -1 width from the near and far 
sides of the circumstellar shell is not resolved with a (J.'Q3 beam 
(Norris et al.Q984). Because the diameter of the OH shell is 
3'.'8 (Bowers & Johnston 1990), the size of the compact area 
unresolved with the VLBI is less than 0.8% of the shell size. 
Therefore it is very likely that the emission from both extreme 
sides of the shell is not depolarized if observed with the NRT 
and it reliably probes a local magnetic field. Small differences in 
the polarization angles at the near and far edges of the OH shells 
in our subsample suggest a regular magnetic field geometry in 
the 1612 MHz maser regions. However, it is hard to explain this 
homogeneous magnetic field over OH shells of large diameters 
(e.g. Habing 1996 ). A plausible explanation for the low values of 
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Fig. 5. Diagram of the circular (mc) and linear (tol) fractional polarizations for the whole sample. The fractional polarization mc 
and hjl are determined for the Stokes V and p peaks, respectively. The sources with a polarized flux density lower than 3<x detection 
level are assumed to have a fractional polarization equal to zero. The points of \mc\ > 20% or/and > 10% are labelled with the 
names of the sources. The thick dashed line defines the relation between the fractional polarizations of the form mc = l.ltn^ + 1.5. 



A^bs is that the origin of the field is intrinsic and the same field 
orientation is carried away (frozen in) by the stellar wind flow- 
ing towards the front and back of the envelope. Alternatively, the 
magnetic field could be of galactic origin and be amplified and/or 
distorted by the stellar wind. High values of A^bs in a minority 
of sources suggest a deviation from globally ordered magnetic 
field caused by local outflows or outbursts. The latter possibil- 
ity is likely for OH17.7-2.0 where Axbr is only -3.7 (Tab.|7) 
but the polarization angle of the eruptive feature near 73 km s -1 
(Szymczak & Gerard2005) is about -15° and differs by about 
64° from the mean value for the extreme blue- and red-shifted 
velocities (Tab.|7]i. 

Nine sources in the sample have linearly polarized emission 
at the two frequencies from the same usually red-shifted side of 
the shell (Tab. A.l, Fig. B.l). The scatters of^- angles within the 
OH mainline channels are on average 2.7 times larger than those 
within 1612 MHz lines. This may suggest that the 1667 MHz 
masers come from more turbulent regions and/or the magnetic 
fields are less ordered. For seven sources the velocities of the 
extreme emission at both frequencies overlap within less than 
0.1 km s" 1 . In the two objects OH127.8-0.0 and OH53.6-0.2 
the difference in the average polarization angle for four neigh- 
bouring channels of the 1612 and 1667 MHz lines is less than 
18° while for the remaining sources it ranges from 36 to 80°. 
The mean difference for nine sources is 44+8° and the median 



is 47°. This result suggests that the magnetic fields probed by the 
two lines are generally not aligned. 

3.6. Possible Zeeman pairs 

The main line V Stokes spectra of RCnc and WHya contain S- 
shaped features possibly caused by the Zeeman effect (Figs. [7] 
and[8]l. Although this assumption needs to be verified by VLBI 
observations one, can roughly estimate the strength of the mag- 
netic field assuming that the ratio of the Zeeman splitting to the 
Doppler linewidth is higher or equal to unity. In this regime the 
separation between the peaks of V profile is a direct measure 
of the field strength (Elitzur 1996). We used the Zeeman split- 
ting coefficient of 0.5886km s -1 mG l and 0.3524km s 1 mG" 1 
for the 1665 and 1667 MHz lines, respectively. For the 1665 
and 1667MHz features of RCnc near 15.2 km s~' the field 
strength estimates along the line of sight are -1.31+0.07mG 
and -2.26+0. llmG, respectively. For the strongest feature of 
WHya near 35.6 km s~' the field strength is -0.32+0.03 mG and 
-0.36+0.04mG at the 1665 and 1667MHz lines, respectively. A 
negative sign means a field pointed towards the observer. These 
estimates should be treated with caution because the Zeeman 
patterns are far from perfect. Nevertheless, in both sources the 
field estimates are quite similar at both frequencies and more- 
over they are consistent with the high angular resolution mea- 
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Table 7. Averaged polarization position angles^ and^R for the 
extreme blue- and red-shifted parts of the 1612 MHz spectrum, 
respectively. In most cases the number of averaged channels is 
four. The A^bk is a difference between^ and^R. The standard 
deviations (SD) are given in the parentheses. 



Name 


Xb (SD) 


Xr (SD) 


AXbr (SD) 


WXPsc 


52.7(4.8) 


63.4(2.6) 


-10.8(5.4) 


OH127.8+0.0 


-20.5(4.4) 


-33.6(3-3) 


13.1(5.4) 


OH138.0+7.2 


-47.3(2.5) 


-34.2(2.3) 


-13.1(3.4) 


IKTau 


15.8(2.1) 3 


30.7(5.2) 3 


-14.8(5.6) 


OH345.0+15.7 


-65.0(1.3) 


-56.1(0.7) 


-8.8(1.5) 


IRAS 17168-3736 


66.6(1.8) 


68.2(1.7) 


-1.6(2.5) 


IRAS 17 177-3627 


32.7(1.6) 


-9.2(1.7) 


41.8(2.3) 


IRAS 17271-3425 


80.6(3.2) 


74.1(2.7) 


6.5(4.2) 


OH1.2+1.3 


79.3(1.0) 3 


80.3(2.5) 


-1.0(2.7) 


IRAS1741 1-3154 


-13.7(0.7) 


-18.7(0.8) 


4.9(1.1) 


OH1. 1-0.8 


-59.8(2.3) 


-83.9(6.4) 


24.1(6.8) 


IRAS17579-3121 


-74.2(6.3) 


-66.1(8.1) 


-8.1(10.3) 


OH15.4+1.9 


-74.3(4.3) 


-89.3(3.9) 


15.0(5.8) 


OH15.7+0.8 


-4.0(1.5) 


-3.1(2.0) 


-0.8(2.5) 


OH12.8-1.9 


-19.5(2.3) 


-9.4(1.1) 


-10.2(2.5) 


OH16.1-0.3 


76.5(3.5) 


90.6(2.4) 


-14.1(4.2) 


OH17.7-2.0 


-76.8(6.0) 3 


-80.4(7.8) 


-3.7(9.8) 


OH25.1-0.3 


44.0(4.8) 


72.6(6.7) 


-28.6(8.2) 


IRAS 18445-0238 


-54.3(4.0) 


-46.2(2.7) 


-8.1(4.8) 


OH32.8-0.3 


15.6(3.2) 


-15.4(10.3) 3 


31.0(10.8) 


OH39.7+1.5 


-57.5(3.3) 


-73.0(3.4) 3 


15.5(4.7) 


OH37. 1-0.8 


-2.3(3.1) 


-13.9(19.2) 2 


11.6(19.4) 


IRAS 19059-2219 


-47.3(3.1) 


-44.2(7.3) 


-3.1(7.9) 


OH42.3-0.1 


38.8(21. 9) 3 


-58.3(4.1) 


-82.9(22.3) 


OH45.5+0.0 


-32.4(7.8) 


17.4(10.1) 


-49.8(12.8) 


OH55.0+0.7 


24.8(1.4) 


20.9(5.4) 


3.9(5.6) 


OH77.9+0.2 


-4.0(4.0) 


67.5(7.5) 3 


-71.4(8.5) 


OH75.3-1.8 


2.4(2.6) 


13.1(1.6) 


-10.7(3.0) 


OH104.9+2.4 


107.7(7.6) 


78.0(2.6) 


29.7(8.0) 


3 three, 2 two averag 
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surements of other AGB and post-AGB stars (e.g. Szymczak et 
al. [1998l Bains et al. l2003l|20U9] l. The presence of polarized fea- 
tures in lots of our sources suggests that many Zeeman pairs will 
be found in high angular resolution studies. 



3.7. OH 1667 MHz overshoot 

There are 63 objects in the sample for which the spectral extent 
of 1667 MHz emission is higher by at least 0.2kms _I than that 
of 1 6 1 2 MHz emission. This overshoot effect is seen either at the 
blue- or red-shifted or at both edges. A summary of the mean and 
median expansion velocities and overshoots for the three classes 
of objects is given in Table [8] Here, the expansion velocity is 
determined as a half of the peak-to-peak velocity width of the 
1612MHz profile. The class of Mira+SR stars has expansion 
velocities significantly lower than in the two other classes. The 
lowest mean and median overshoots of 0.78±0.16kms _1 and 
0.47 km s~', respectively, belong to the Miras and SR objects. 
The highest median value of overshoot of 2.08 km s ) is deter- 
mined for the post-AGB+PN class. The phenomenon is most 
frequent (63%) in the OH/IR objects. 

The overshoot effect was broadly discussed by Sivagnanam 
& David ( fT999l > for a sample of 21 AGB and post-A GB ob- 
jects who proposed several mechanisms. Deacon et al. (2004) 
suggested the existence of some acceleration in the outer parts 
of the shell where the 1612 and 1667 MHz maser originate and 
that the 1667 MHz masers are located farther than the 1612 MHz 




Velocity (km s" ) 



Fig. 6. Examples of the 1612 MHz p spectra in the sources with 
the emission appeared at the blue- and red-shifted velocities. The 
polarization position angles (points) are superimposed (right or- 
dinate) for the extreme outside channels. The velocity is mea- 
sured with regard to the systemic velocity. 

ones. The trend in overshoot found in our sample suggests that 
the phenomenon can be related to the evolutionary sequence. We 
did not find any obvious relationships between overshoot and the 
polarization parameters. 

4. Discussion 

4.1. Occurence of polarization 

The present full-polarization data for a substantial sample of 
AGB and post-AGB objects have allowed us to give a com- 
prehensive picture of the polarization properties of OH circum- 
stellar masers. First of all, our high sensitivity and spectral res- 
olution study shows that the polarized OH maser emission is 
widespread. This finding is at odds with early reports, which 
demonstrated a paucity of polarized features in the 1612 MHz 
spectra and only occasional detections of linear polarization 
at 1667 MHz (Cohen [T989l for review). Thus, spectral resolu- 
tion appears to be a crucial factor (Cohen et al. 119871 Zell & 
Fix fTWTT l. 

The study allowed us to extract a complete subsample of 
1612 MHz emitters, composed of 8 Miras+SR stars, 38 OH/IR 
objects and 10 post-AGB+PN sources, which is not biased by 
the sensitivity limit and the minimum believable fractional po- 
larization (see Sect. 3.2 for definition). In this subsample the de- 
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Table 8. Average and median values of expansion velocities and overshoots for the three classes of sources. 



Class 


Line 
(MHz) 


Number 


Expansion velocity (km s 1 ) 


Number 


Overshoot (km s ') 


Mean(SD) 


Median 


Mean(SD) 


Median 


Mira+SR 


1612 


23 


10.68(1.08) 


12.16 










1667 


26 


8.36(0.98) 


6.12 


16 


0.78(0.16) 


0.47 


OH/IR 


1612 


56 


15.15(0.57) 


14.38 










1667 


37 


15.36(0.54) 


15.09 


37 


1.18(0.14) 


0.93 


post-AGB +PN 


1612 


14 


12.08(1.14) 


12.36 










1667 


8 


13.34(0.74) 


13.51 


10 


2.25(0.62) 


2.08 
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Fig. 7. Tentative Zeeman splitting in the OH maser lines of 
RCnc. The left- (LHC) and right-hand (RHC) circular polar- 
ization spectra, V Stokes and p spectra are shown for the two 
frequencies. 
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Fig. 8. Same as in Fig.|7]but for W Hya. 



tection rate of polarized features is 75-90% and there is no intrin- 
sic difference in the occurrence of polarized emission between 
the three classes of studied objects. A complete subsample of 
1667 MHz sources contains only 2 Miras+SRs, 3 OH/IR objects 
and 4 post-AGB+PN stars. All but one of these targets show po- 
larized emission. It seems that the occurrence of polarized fea- 
tures at 1667 MHz does not differ between the three classes of 
stars. 

4.2. Effect of depolarization 

For 1612 MHz sources the mean and median degrees of circu- 
lar polarization in the post-AGB+PN objects are significantly 



higher that those in the two other classes of studied stars. We 
notice that in the sample are at least two bursting post-AGB 
stars; IRAS16342-3814 and OH17.7-2.0 with m c >80%. These 
objects sometimes experience OH bursts of highly polarized 
emission at 1612 MHz usually not associated with bursts at 
1665 and 1667 MHz (Szymczak & Gerard 2005). In many post- 
AGB objects there is a disruption of detached AGB 1612 MHz 
shells that are no longer spherically symmetric and the appear- 
ance of a well-defined bipolar morphology at 1667 MHz is ob- 
served (Zijlstra et al. 120011 . In these partially-filled 1612-MHz 
structures, observed at high angular resolution, Zeeman com- 
ponents are much less spatial blended than in dense (10~ 4 Mq, 
Habing 1996), more fully filled shells that are observed around 
OH/IR objects (e.g. Bowers & Johnston 1990). Therefore, the 
enhancement of the circular polarization degree at 1612 MHz in 
the post-AGB+PN objects compared to the two other classes of 
objects could be caused by negligible depolarization in the first 
class of objects. However, the post-AGB+PN group may be het- 
erogeneous because of the lack of accurate criteria for selection 
(Sect. 2.1) and significant depolarization cannot be excluded in 
some objects (Szymczak & Gerard 2004). 

The mean degree of circular polarization at 1667 MHz is sig- 
nificantly higher in the Mira+SR stars than in the other classes 
of objects. The mass loss rates for Mira and SR variables are 
about 10~ 7 M Q (Habing [1996l ) while the 1667 MHz maser orig- 
inates in the inner parts of the shell at a distance of about 50- 
150 AU from the central star (Chapman et al. l 19941 Szymczak 
et al. l 19981 ). This implies that regardless of the magnetic field 
type, its strength in the 1667 MHz maser regions should be 
higher than in the 1612MHz OH maser outer shells. High an- 
gular resolution images revealed a highly fragmentary distribu- 
tion of the 1667 MHz emission (Chapman et al. l 19941 Szymczak 
et al. l 19981 ). indicating a low filling factor for the OH molecule. 
In this shell the spatial overlapping of Zeeman components is 
weak, if there is any at all. Furthermore, Miras and SR variables 
often exhibit erratic variability of OH emission, which usually is 
highly polarized (Etoka & Le Squeren 1997). We conclude that 
the higher degree of circular polarization in the Mira+SR stars 
compared to the other classes of objects is because of negligible 
depolarization in the envelopes that are partially filled by the gas 
clouds sustaining OH maser emission. 

4.3. Zeeman splitting 

The interpretation of OH maser polarization data strongly de- 
pends on the assumed value of the ratio of Zeeman splitting to 
Doppler line width, xb, (Elitzur 1996). For xb > 1 the separa- 
tion of the cr components of the Zeeman pair scales directly 
to the strength of the magnetic field along the line of sight. If 
xb < 1 only the Doppler width can be estimated from the sep- 
aration of the two peaks of V Stokes profile. High angular res- 
olution observations of OH lines in some AGB and post-AGB 
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stars (Szym czak et al. [T998l Bains et al. l2003ll2004ll2009l Amiri 
et al. l20101 > revealed a magnetic field strength of 1 -3 mG. This 
value agrees well with our tentative estimates of magnetic fields 
of 0.3-2.3 mG in two Miras. The kinematic temperature in re- 
gions of circumstellar 1612 and 1667 MHz masers is 50-100K 
(Goldreich & Scoville [T976l Fong et al. 120021 hence for OH 
molecule the thermal linewidth is 0.36-0.51 km s . It implies 
that > 1 is fulfilled for most sources in the sample. We note 
that the Zeeman pairs are of unequal strength, and in several 
cases likely only one sense of circularly polarized emission is 
seen. This can occur when the maser radiation is amplified in 
the region where the gradients of the velocity match the mag- 
netic field (Cook 1966 ). A simpler explanation for uneven circu- 
lar polarization is the overlap of different Zeeman components 
caused by a velocity gradient within the maser region (Deguchi 
& Watson [l986l 

In the case of xr > 1 the polarization of the <x com- 
ponents is described by the equations Vy7=2cos0/(l+cos 2 f9) 
and e//=sin 2 <9/(l+cos 2 0) (Goldreich et al. [l973l Elitzur [T996l 
where 8 is the angle between the magnetic field direction and 
the line of sight. We used these equations to model the ;«c - 
diagram. It is easy to show that the ratio V/I integrated and aver- 
aged over the solid angle builds up rapidly between = 0° and 
50° while the Q/I has barely started growing. This means that 
the ratio V/Q is very high until 45° and that may explain the 
void in Fig. [5]even with a spherically symmetric distribution of 
the 8 angle. Although our study was of insufficient angular res- 
olution to determine polarization properties of individual maser 
spots and the strongest polarized features were used to construct 
the ntc - m-L diagram, the appearance of the void (Fig. [5} may be 
the proof that the polarized features are the cr components and 
that the condition x% > 1 is fulfilled for unsaturated OH masers. 

Fig.[5]shows a distinct group of 1612/1667 Mira+SR sources 
with only circularly polarized emission (the first six objects in 
Tab. indicating that the magnetic field is parallel to the line 
of sight. Another small group of mainly OH/IR objects (the five 
last objects in Tab. [6} with predominant linear polarization at 
1612 or 1667 MHz represents the case where the magnetic field 
is perpendicular to the line of sight. These objects appear to be 
good candidates for high angular resolution studies to determine 
their properties better and to examine if significant differences 
in polarization parameters could be caused by different states of 
saturation as predicted in the standard model of maser polariza- 
tion (Goldreich et al. [1973l Vlemmings l2006t . 

5. Conclusions 

We have presented the results of full-polarization observations of 
the 1612 and 1667 MHz OH maser lines towards a sample of 1 17 
AGB and post-AGB stars. A complete set of full-polarization 
spectra and the basic polarization parameters was given. 

The polarized features occur in more than 75% of the sources 
in the sample complete to the total flux density higher than 7.5 Jy 
and the minimum fractional polarization of 2%. There is no in- 
trinsic difference in the occurrence of polarized emission be- 
tween the three classes of objects. The highest fractional polar- 
ization occurs for the post-AGB+PN and the Mira+SR classes 
at 1612 and 1667 MHz, respectively. Differences in the frac- 
tional polarizations are likely caused by depolarization caused 
by blending of polarized emission. 

The relationship between the circular and linear fractional 
polarizations for extreme red- and blue-shifted features, which 
are less likely to be blended, was interpreted within the frame of 
the standard models of polarization for the regime of the Zeeman 



splitting higher than the Doppler line width. The strength of the 
magnetic field derived from the tentative Zeeman pairs for a few 
Miras of 0.3-2.3 mG is consistent with the results of high angular 
resolution observations of other AGB stars. The front and back 
caps show only small differences in polarization angle (<20°) 
at 1612 MHz, which suggests a regular geometry of the mag- 
netic field in the outer parts of the shells. The observed relation 
between the circular and linear fractional polarizations suggests 
that the polarized features are the cr components and the Zeeman 
splitting is higher than the Doppler line width. 

The study has provided many targets for future investiga- 
tions. VLBI + e-MERLIN mapping of polarized OH maser emis- 
sion will be essential to measure the magnetic field strength and 
its structure and to refine the present results. 
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